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Abstract 

We analyze the microscopic, topological structure of the interface between domains 
of opposite magnetization in 3D Ising model near the critical point. This interface 
exhibits a fractal behaviour with a high density of handles. The mean area is an 
almost linear function of the genus. The entropy exponent is affected by strong finite- 
size effects. 
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1. Introduction 

The interface in 3D statistical systems above the roughening temperature behaves like a 
free, fluctuating surface with a topology determined by the boundary conditions. Recently, 
new computational algorithms have been applied to these systems [1-6], providing us with 
a powerful tool for testing our ideas on the behaviour of fluctuating surfaces. Alternatively, 
these surfaces can be thought of as space-time histories of closed strings which can be used 
to describe the infrared behaviour of the dual Z2 gauge theory. 

In a previous paper |?J we studied the free-energy of the interface as a function of 
its shape. We found rather strong shape effects which are accurately described by the 
gaussian limit of the Nambu-Goto action. 

This observation suggests that interface configurations are mostly made of smooth 
surfaces subjected to long- wavelength fluctuations which account for the observed finite- 
size effects. On the other hand, this picture seems to strongly disagree with the process 
of crumpling, which should take place for random surfaces embedded in a target space of 
dimension d > 1. In other words, if the interface is described by the Nambu-Goto action, 
it should take the shape of a branched polymer, rather than that of a smooth surface. 

In this letter we face this dilemma by studying the interface from a microscopic point 
of view. Actually, it turns out that the interface at small scales is much more similar to a 
sponge than to a smooth surface. In particular there is a strong, almost linear correlation 
between the area and the genus of the surface, indicating the formation of a large number of 
microscopic handles which is proportional to its area. Yet, the partition function summed 
over all the genera behaves like that of a smooth toroidal surface: this means that the huge 
number of microscopic handles, produced by the instability toward crumpling, have as the 
net effect a simple non-perturbative renormalization of the physical quantities associated 
to the surface, according to an old conjecture on string theory [||]. 



1 



2. The method 

The first problem to be solved in order to study the microscopic structure of the 
interface is to express the 3D Ising model or its dual gauge version as a gas of suitable 
closed surfaces. There are essentially two different ways to do it. One is based on the 
strong coupling expansion of the dual gauge model |]] and involves also self-intersecting 
and non-orientable surfaces. Here instead we identify the surfaces as the Peierls interfaces 
of an arbitrary Ising spin configuration on a cubic lattice By construction these are 
closed, orientable surfaces composed of plaquettes of the dual lattice orthogonal to each 
frustrated link. In this way each plaquette appears at most once in the construction of 
the surface S. 

Tab. I 

Graphs for the microscopic reconstruction of the interface. Each cube is dual to a vertex V q with 
q edges of the interface as drawn in fig. 1. The encircled sites have a spin different from the other 
sites. Graphs with different vertex decompositions correspond to end-points of contact lines. The 
first decomposition is for positive vertices, while the other is for negative vertices. 
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A link belonging to S is said to be regular if it glues two plaquettes of S. The only- 
singularities which may appear on these surfaces are links glueing four distinct plaquettes 
of S. These singular links form contact lines of the surface which are sometimes improperly 
called self-intersections. These contact lines are the only potential sources of ambiguities 
in the topological reconstruction of the surface. We shall see shortly how to remove them 
in a simple, consistent way. A vertex of the surface S is dual to an elementary cube of 
the lattice of the Ising configurations. According to the distribution of spins inside the 
cube, there are twelve distinct vertices as listed in Table I. Some of them include singular 
links and may be interpreted as the coalescence of distinct, regular vertices. Only the end- 
points of contact lines can be decomposed into two or more inequivalent ways: for instance 
the first graph in the second row of the table can be thought either as the coalescence of 
two vertices with three edges, or as a single vertex with six edges and the singular link is 
splitted accordingly into two different ways as drawn in figure 1. 




Fig. 1 

Example of decodification of a graph of Table I. This cube has six frustrated links, labelled by 
a, b, . . . , /, which are dual to six plaquettes forming a vertex with a singular link. It can be 
considered either as a 6-vertex, or as two 3-vertices. This choice depends on the sign of the 
magnetization of the cube, which can be in this case ±4. 

Choosing arbitrarily this splitting procedure for each end-point may generate global 
obstructions in the surface reconstruction, due to a constraint which links together the 
splitting of two end-points connected by a contact line. In order to formulate this con- 
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straint, it is useful to attribute a sign to each singular vertex according to the sign of the 
magnetization of the corresponding elementary cube. Then, working out some explicit 
example, it is easy to convince oneself that there are no global obstructions if one chooses 
different splittings for vertices of different sign. More precisely, the only constraint to be 
fulfilled in the replacement of each singular link with a pair of regular ones is that two 
end-points connected by a contact line must be splitted in the same way if they have 
the same sign, or in the opposite way if their sign is different. The two decompositions 
appearing for singular vertices in table I correspond precisely to the two possible signs 
of the vertex. The only case in which the surface reconstruction is not immediate is the 
first graph of the last row of tab. I, corresponding to a vertex with six singular lines: the 
three possible decompositions cannot be selected by the sign of its magnetization (which 
is zero), but by the signs of the end-points of the six contact lines. 

The main consequence of the above construction is that we have obtained a simple rule 
to assign unambiguously to each Ising configuration a set of self-avoiding closed random 
surfaces. Notice that this fact agrees with a result of David [|ll]] who found that a gas of 
self-avoiding surfaces in a special three-dimensional lattice belongs to the same universality 
class of the Ising model. 

We may now evaluate the genus h of each self-avoiding surface through the Euler 
relation 

F-E + V = X {S) = 2-2h (1) 

where F is the number of faces (plaquettes), E the number of edges, V the number of 
vertices and x t ne Euler characteristic; the genus h gives the number of handles. Denoting 
by N q the number of vertices of coordination number q, we have obviously, according to 
table I, 

V = N 3 + N± + N 5 + N 6 + N 7 ; 4F = 2E = 3N 3 + 4iV 4 + 5N 5 + 6iV 6 + 7N 7 , (2) 
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which gives at once an even simpler expression for the genus of this kind of surfaces 

4(y -F) = N 3 -N 5 - 2N 6 -3N 7 = 8-8h . (3) 



3. Results 

We are interested on the topology of the interface between two macroscopic domains of 
opposite magnetization. For this reason we consider very elongated lattices with periodic 
boundary conditions in the two short directions (denoted in the following by L) and 
antiperiodic boundary conditions in the long direction (denoted by L z ). This forces the 
formation of an odd number of interfaces in the L z direction. We then isolate one of these 
interfaces by reconstructing all the spin clusters of the configuration, keeping the largest 
one and flipping the others. Note, as a side remark, that for the values of (3, L and L z 
we have studied, typical configurations contained only one macroscopic cluster, besides 
a huge number of microscopic ones. It is now easy to evaluate area of the interface by 
counting simply the number of frustrated links of this cleaned configuration. 

The growth of this area as a function of the size of the lattice gives a first description 
of the fractal behaviour of the interface: in a set of Monte Carlo simulations near the 
critical point on elongated lattices L 2 x L z with L z > 120 and 8 < L < 16 (see table II) 
we found that the mean area < F > of the interface is a strongly varying function of the 
transverse size L, well parametrized by the following power law 

<F>=nL dH , (4) 

with du ~ 3.7 and k ~ 0.47 . As a consequence, a sizeable fraction of the lattice is invaded 
by the interface, and a small increasing of the transverse lattice section implies a rapid 
growth of the area of the interface, so we have to take very elongated lattices in order to 
avoid wrapping of the interface around the antiperiodic direction, which would give rise 
to unwanted finite volume effects. 
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Each interface S of area F of an arbitrary Ising configuration contributes to the par- 
tition function simply with a term e~^ F ; then we can define the following generating 
functional Z 

m = E E e_/3F = EE z h( F ) > ( 5 ) 

F # surfaces of area F F h 

where Zh(F) is the partition function for a surface of area F and genus h. Since we are 
dealing with macroscopic surfaces, we can consistently assume that their multiplicity is 
well described by the asymptotic behaviour of the entropy of large random surfaces fll 
which yields 

ZhiF)^^ ~ F 6 *- 1 e^ F e- /JF = F^V^ , (6) 

where, using the terminology of two-dimensional quantum gravity (2DQG), fi is the cos- 
mological constant and the exponent b x is a function of the Euler characteristic \- I n the 
continuum theory of 2DQG this exponent can be evaluated exactly for the coupling to 
conformal matter of central charge c < 1 [p^| . The result is 

b x = - b - X = b{h-l) , (7) 

with 

= 25-c+ V(l-c)(25-c) 

12 1 ' 

b is a monotonically decreasing function of c for c < 1 and becomes complex for c > 1, 

where the surface get crumpled and these formulas lose any physical meaning. It is known 

that for c = 1 there are logarithmic corrections to eq.(||), while for c > 1 it is only known 

that the number of surfaces of given area is exponentially bounded [12|] . It turns out that 

eq.(||) fits well to our numerical data. 

We can evaluate the exponent b x for the interface by measuring the mean area at fixed 

genus. Indeed defining 

Z h { l x)=Y J Z h (F) , (9) 
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we have 

<F> h =—?-logZ h (fi)~^ . (10) 

In figure 2 we report a typical outcome of this analysis. Actually < F >h is dominated 
by a linear term, like in eq.(R), but there is also a small, negative quadratic contribution 
which probably accounts for the self-avoidness constraint: when the number of handles 
h is very large the excluded volume effects (proportional to h 2 ) become important. The 
data are well fitted by 

^ a b , d , 9 , . 

< F > h = — | — h H — h . (11) 
H n fx 

The slope b can be called the entropy exponent, and a is related to the strings susceptibility 
by a = 7 S — 2. Comparing eq.([Tll) with eq.([|) we have of course, for c < 1, d = and 

b = -a = 2 - 7 S . (12) 

Since we are dealing with surfaces of very high genus h ~ 10 2 , we can evaluate bj\i 
and d//j, ~ — 10~ 4 b//j, very accurately, while the constant a is affected by large errors and 



cannot used to evaluate j s , nevertheless we shall argue shortly that for our surfaces eq.(12) 
is no longer true. 



Combining eq.(10) with other physical quantities at fixed genus it is possible to evaluate 
the cosmological constant [i as a function of (3 and of the lattice size. In particular, using 
the derivative dZ gjf^ , we get easily the following equation 

log(l//i) = -V>(6 x )+<logF> , (13) 



where ip denotes the logarithmic derivative of the Euler T function. This formula fits well 
the numerical data for large areas and we used it to evaluate /i. The results are reported 
in table II as well as the entropy exponent b defined in eq.([TT|) . It turns out that b 
is affected by rather strong finite size effects, being a decreasing function of the lattice 
size L an hence of the area. Similar finite size effects have been observed for the string 
susceptibility exponent in the planar quantum gravity coupled to c > 1 matter [0] . Note 



however that, if one assumes eg. (|12|) , one should conclude that 7 S exceeds the theoretical 
upper bound 7 S = 



15(1 . A possible way out is that eq. (|12|) is not fulfilled by this kind 



of surfaces, indeed the argument leading to such upper bound assumes h = and gives 



no restriction on b. On the other hand there is a renormalization group argument [16 
showing clearly that eq. (|l2|) is justified only for c < 1. 



Tab. II 

Mean area < F >, entropy exponent b, cosmological constant n and renormalized 
cosmological constant fiji of the interface in a set of elongated lattices of shape L 2 x L z at 
p- j3 c = 0.0059 . 



L 2 x L z 


< F > 


b 






8 2 x 120 


1057(10) 


0.73(20) 


0.0141(39) 


0.00122(4) 


10 2 x 120 


2434(19) 


0.61(15) 


0.0097(24) 


0.000613(15) 


12 2 x 180 


4719(47) 


0.30(3) 


0.0047(5) 


0.000347(4) 


14 2 x 240 


8506(103) 


0.19(3) 


0.0026(5) 


0.000208(3) 


16 2 x 320 


14803(249) 


0.11(2) 


0.0015(2) 


0.000125(5) 



It is also possible to define a renormalized cosmological constant ^lr by considering 
the sum over all the genera Z(F) = J2h^h(F) . It turns out that for large areas the 
distribution of surfaces is accurately described by an exponential fall off of the type 



Z(F) F ^ 00 ~ ce 



(14) 



as shown in fig. 3. Comparing this equation with eq.(^), we argue that the sum over all the 
topologies produces as net effect an effective interface which behaves as a smooth surface 
with a different cosmological constant. Actually the renormalization effect is very large: 

is about one order of magnitude smaller than the unrenormalized quantity fj, (see 
table II). This phenomenon seems the two-dimensional analogue of a quantum gravity 
effect described by Coleman |l7|], who argued that the sum over topologies has the effect 
of making the cosmological constant vanishing small. 

In conclusion, we have found a simple way to generate high-genus self-avoiding random 



surfaces. The most interesting question is of course whether their study will lead to new 
insights in the physics of random surfaces. 
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Figure Captions 

Figure 2 The mean area of the interface as a function of the genus in a lattice of size 
12 2 x 240 at - (5 C = 0.0044 . 

Figure 3 The multiplicity of interface configurations summed over all the topologies as a 
function of the area in a lattice of size 14 2 x 240 at (3- (3 C = 0.0059. The lower set 
of data is the multipicity a fixed genus h = 70. The stright line is the fit to eq . (|i~4|) . 
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